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Doping dependence of the exchange energies in bilayer manganites: Role of orbital 

degrees of freedom 
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Recently, an intriguing doping dependence of the exchange energies in the bilayer manganites 
La2-2xSri+2iMn207 has been observed in the neutron scattering experiments. The intra-layer 
exchange only weakly changed with doping while the inter-layer one drastically decreased. Here 
we propose a theory which accounts for these experimental findings. We argue, that the observed 
striking doping dependence of the exchange energies can be attributed to the evaluation of the 
orbital level splitting with doping. The latter is handled by the interplay between Jahn- Teller effect 
(supporting an axial orbital) and the orbital anisotropy of the electronic band in the bilayer structure 
(promoting an in-plane orbital), which is monitored by the Coulomb repulsion. The presented theory, 
while being a mean-field type, describes well the experimental data and also gives the estimates of 
the several interesting energy scales involved in the problem. 
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The three-dimensional (3D) cubic manganites 
Ri-^B^MnOa (where R is trivalent rare-earth and B 
is divalent alkaline ion, respectively) as well as the two- 
dimensional (2D) bilayer ones La2-2£cSri + 2a;Mn207 have 
attracted recent interest not only due to the discovery of 
colossal magnetoresistance (CMR) in these compounds 
but also because of their rich and rather unusual phys- 
ical properties B At different doping concentration the 
manganese oxides exhibit a wide diversity of the ground 
states, and the small changes of some external param- 
eters of the system such as doping, chemical pressure, 
temperature or magnetic field can result in a drastic 
modification of the physical properties as well as cause 
the transition from one to another ground state. 

In fact, in very recent spin- wave measurements on 
the bilayer manganites La2_2a;Sri+2a;Mn207 (in the dop- 
ing range x = 0.3 — 0.5) strongly anisotropic doping 
dependence of the inter-layer (within the bilayer) (Jj_) 
and of thc_intra-layer ( Jii ) exchange couplings has been 
revealedm! At x = 0.3 both exchange constants are fer- 
romagnetic,-(FM) and of the same order J±=5 meV, and 
J||=4 meVE With further doping (x > 0.3), Jm changes 
very slightly while inter-layer exchange rapidly decreases, 
by a factor of four at x = 0.4, and changes sign at 
x ~ 0.45. 

Within the double-exchange (DE) picture FM cou- 
pling between the nearest-neighbour (NN) tig core spins 
is mediated by the hopping of e 9 -electrons (d x 2_ y 2 and 
d 3z 2_ r 2 orbital states) and scales as a kinetic energy of 
these electrons on a given bond. Since only d 3z 2_ r 2 (ax- 
ial) orbital has a finite inter-layer transfer amplitude (z- 
axes is set perpendicular to a bilayer), the inter-layer 
exchange is only meditated by this orbital. While the 
intra-layer one is mainly determined by the d x 2_ y 2 (in- 
plane) orbital state, that has a highest intra-plane trans- 
fer amplitude. Therefore, above discussed decrease of 



the inter-layer exchange constant with doping can be 
ascribed to the change of the nature of occupied or- 
bital state from mainly d 3z 2_ r 2 character, at low d 
ing x < 0.3, to mostly d x 2_ y 2 one at higher doping! 
This qualitative picture is also consistent with the recent 
x-ray magnetic Compton-profile measurements. □ The re- 
sults of this experiment have shown that upon doping 
(for 0.35 < x < 0.45 ) the electrons are mainly with- 
drawn from d3 Z 2_ r 2 orbital state and the occupancy of 
d x 2_ y 2 state remains practically the same. However, we 
point out, that the observed doping dependence of the ex- 
change constants and of the orbital occupancy, being in- 
terrelated, can not be understood within a rigid band pic- 
ture that assumes a constant value of orbital level split- 
ting. Interpretation of the experimental findings within 
this picture would imply that only d^ z 2_ r 2 band crosses 
the Fermi level and is emptied out upon doping, while 



the ef. 



band lies below and shows a constant filling. 



From the other hand, if d x 2_ y 2 band does not cross Fermi 
surface, then it is no longer active to generate the intra- 
layer ferromagnetic DE interaction leading to Jii < J±, 
in contrast to experimental findings Jn > J±. 

The aim of the present work is to propose a mecha- 
nism that can explain the experimental findings. The 
paper is organized as follows: First we briefly discuss 
the structural aspects of the compound relevant for our 
study. Next we give the basic physical ideas and set 
up the minimal model Hamiltonian to provide a micro- 
scopic description. We then present our results on the 
doping evaluation of the orbital occupancies and the ex- 
change energies, and finally, compare the theoretical find- 
ings with that of the experiments. 

The compound La2-2a;Sri + 2xMn207 consists of bi- 
layer slices of MnOg octahedra, separated by insulating 
(La,Sr) 2 02 layers that serve to decouple the bilayers both 
electronically and magnetically. Therefore, as a first ap- 
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proximation the system can be treated as composed of 
independent bilayers. One electron placed in the doubly 
degenerate e g level, makes the Mn +3 0^ 2 complex Jahn- 
Teller ( JT) active. The experimental analysis of the crys- 
tal strutter has revealed, that JT distortion results in the 
elongation of MnC>6 octahedra along z-axis.U Five bonds 
(four equatorial and one apical with shared oxygen ion) 
are essentially identical. The elongation of octahedra oc- 
curs mainly because the apical bond with unshared oxy- 
gen atom is longer than the other five bonds. These type 
of JT distortion promotes the axially directed d 3z 2_,,2 or- 
bital by lowering its on-site energy. The dimensionless 
parameter describing the static Jahn- Teller distortion 



can be defined as <5 JT = ((dX-o> - Wn-oll/iC'' 
where (d^}^ 
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is an averaged apical (equatorial) Mn- 
O length. The relative distortion 5jt also gives the scale 
of JT induced orbital level splitting. 

The physical picture we shall be based on is as fol- 
lows: Both the orbital occupancy and exchange ener- 
gies strongly depend on the orbital level splitting. The 
latter consist of the two parts: JT induced gap and 
correlation induced one. The JT distortion is largest 
at low doping and monotonically relaxes upon doping 
((5 JT ~ 0.036-0.020, and 0.0 for x = 0.3, 0.4, and 0.5, 
respectively) Therefore, at low doping the axial orbital, 
promoted by the JT distortion is predominantly occu- 
pied. The matrix element of the inter-orbital Coulomb 
energy does not depend on the character of predom- 
inantly occupied orbital state. Therefore, on-site re- 
pulsion further supports already preferred orbital state 
and by enhancing the orbital polarization reduces the 
Coulomb energy.Q At low doping, the physics is mainly 
determined locally and kinetic energy term is less impor- 
tant. However, upon increasing the carriers number the 
JT induced gap decreases while the kinetic energy starts 
to play more dominant role. In the bilayer system, the 
latter promotes the in-plane orbital state, that has the 
highest intra-layer transfer amplitude and by forming the 
wider band can lower the kinetic energy of the system. 
Therefore, one would expect that at some doping x = x c 
the possible energy gain due to the kinetic energy will 
overcome the crystal- field splitting of the orbital levels. 
Therefore, for x > x c , it becomes energetically more fa- 
vorable to lower the energy of the in-plane orbital and 
hence, correlation induced splitting will change the sign. 
As will be shown below, the doping dependence of the 
orbital level splitting, resulted from the discussed inter- 
play, gives the consistent description of the experimental 
findings. 

The minimal model Hamiltonian which describes 
above discussed scenario consists of the two contributions 
H = H c \ + H sp , where H e \ and H sp describe charge and 
spin degrees of freedom, respectively.D We retain only the 
preferred spin component of the fermionic operators, in 
the fully polarized (semi-metallic) FM case and also as- 
sume that all the other degrees of freedom are integrated 
out to give effective model parameters relevant to the low 



energy physics. The electronic part of the Hamiltonian 
can be written as: 



dL,d,0 + H.c. 



[(-1) Q A JT - fi] [n ia + ri ia ] U cS \nani2 + nnn l2 ] 

(i) 

The first and the second term of Eq. ([!]) describe an intra- 
and inter-layer electron hopping between the two e g or- 
bitals of the NN Mn-ions, respectively. Index i numbers 
the unit cell composed by two Mn sites, di a and di a are 
the electron annihilation operators on these two sites. 
The orbitals d^ z 2_ r 2 and d x 2_ y 2 correspond to a(f3) = 1 
and 2, respectively. Intra- and inter-layer transfer matrix 
elements are given by 
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The electron density operator on a given orbital state is 
denoted by rii a , [i is chemical potential, Ajt = 3<5jt is 
JT induced orbital level splitting, g is properly normal- 
ized coupling constant of electrons with JT active phonon 
modes, and <5jt is a measure of JT distortion discussed 
above. In the present paper we do not attempt to calcu- 
late JT distortion by minimizing the total energy of the 
system (including lattice elastic energy) and model the 
dimensionless parameter 5jt by the following doping de- 
pendence <5jt = 0.17(0-5 — x), that reasonably reproduces 
the experimental data.El 

The last term in Eq.(|l|) represents an effective, rel- 
evant for the low energy physics, on-site inter-orbital 
Coulomb repulsion between electrons. We assume that 
U e ff has been already renormalized due to the many- 
body effects and properly screened in the metallic state. 
Therefore, one expects that [/ c ff is. much smaller than 
the bear ionic value of Hubbard Uu We treat this term 
within the mean-filed (MF) approximation by introduc- 
ing the MF parameter describing the anisotropy of the 
orbital occupancy Sn = (n±) — (712) ■ This introduces 
additional splitting of local orbital levels due to the 
electron-electron interaction. Therefore the total split- 
ting of orbital levels will be determined by both JT ef- 
fect and electron-electron repulsion with following value: 
2A = 2Ajt + UrffSn. We emphasize, that the nonzero 
value of MF parameter Sn does not imply any symme- 
try breaking long range orbital order and the symmetry 
of Hamiltonian remains tetragonal. Even in the non- 
interacting case and absence of any JT distortion one 
finds Sn ^ 0. The band structure with tetragonal sym- 
metry brakes the local cubic symmetry and promotes pla- 
nar d x 2_ y 2 orbital with wider band leading to Sn < 0. In 
the orbital pseudo-spin language the ground state can be 
represented by nonzero value of the axial component of 
the orbital pseudo-spin (t£) — ((rin — n i2 )/2) 7^ due 
to the nonzero pseudo-magnetic field produced by the 2D 
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structure of the system, and with no order in the pseudo- 



0. 



spin basal plane {r i 

The resulted MF Hamiltonian is bilinear and 
can be diagonalized by two subsequent canonical 
transformations. Q First step is to introduce the bond- 
ing and anti-bonding states as a(6) a (/3)k = [d a (p)k T 
^a(/3)k]/ 1 /2 and decouple the MF Hamiltonian in two 
parts each consisting of the two orbitals. The next trans- 



formation A(i?)ik 
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izes the two-band Hamiltonian bringing it into the form 
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factors and eigen-frequencies are given by 
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with the following notations c^mjj = — ^7k ± t — /i + A, 

4 2 = - 3i 7k - A* — A, e^ 2 = -V3t7 k , and 7 ( 7 ) k = 
[cos k x ± cosfcyJ/2. 

The above introduced MF parameter Sn and chemical 
potential fj, for a given doping can be obtained from the 
system of self-consistent equations which in terms of the 
obtained eigen-states of MF Hamiltonian reads as: 



(4) 
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where n(e) is the Fermi distribution function. 

Lets us now discuss the spin degrees of freedom of 
the system. In the DE exchange limit (Hund's cou- 
pling Jh 3> W carriers band-width) the spin subsystem 
in the two orbital model can be also mapped to an ef- 
fective NN Heisenberg model H = — £\ g, J\\&iSi+5\\ — 
J2i 6 ± J±Si&i+8± with effective intra-layer J|| and inter- 
layer Jj_ FM exchange couplings defined as Jy = J® E — J, 

J± = J EE — J, and J is the antiferromagnetic super- 
exchange constant between the t2 g spins. Ferromagnetic 
intra- and inter-layer DE energies are given by 
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With the above formulated scheme we proceed as fol- 
lows. First, we solve the MF equations (Q) to determine, 
the orbital level splitting A and chemical potential /i for 
a given doping. Then, expressing the Eq. |^ for the ex- 
change constants in terms of the eigen-states of the MF 
Hamiltonian we end up with the intra- and inter-layer 



exchange energies for a given doping. Model parameters 
are fixed in a way to reproduce the experimental results. 

In Fig. |l| the calculated doping dependence of intra- 
and inter-layer exchange couplings is presented together 
with the experimental data from Refs. (see also 

Rcf. 11). The best fit to the data has been achieved 
for t = 0.18 eV, U cS = 0.7 eV, g = 0.5 eV, and 
25 J =11 meV (see also discussion below). The hop- 
ping integral sets the overall energy scale. The value 
of electron-phonon coupling strength fixes the doping at 
which J|| = Jj_ and system shows the isotropic behavior. 
The interaction term J7 e ff determines the steepness of the 
drop of inter-plane exchange constant and the AFM su- 
perexchange J shifts rigidly the whole picture relative to 
the Y axis of Fig. 1. Therefore, despite the number of in- 
dependent parameters all of them can be unambiguously 
extracted from the fitting. 




0.25 0.3 0.35 0.4 0.45 

Hole concentration x 

FIG. 1. Inter- (SJ_i_) and intra-layer (SJ\\) exchange con- 
stants as a function of hole concentration (solid and dashed 
line, respectively). The filled (open) circles, squares, and di- 
amonds are inter-(intra-)layer exchange constants from neu- 
tron scattering data of Refs. 2, 3, and 9, respectively. 

As it is seen from Fig. [I] the intra-layer exchange is 
practically unaffected by doping while the inter-layer one 
is dramatically reduced. When inequality J^ E < J be- 
comes valid, the antiferromagnetic super exchange pre- 
vail ferromagnetic DE signaling the instability of the FM 
ground state of the bilayer. In fact, neutron diffrac- 
tion stuchf on x = 0.5 sample revealed the A- type AFM 
ordering 

As it follows, the above discussed doping dependence 
of the exchange constants is in one to one correspon- 
dence with the doping dependence of the orbital occu- 
pancy. The latter is shown in Fig. 2 calculated for the 
same values of the model parameters. The experimental 
data from Ref. |J is also presented on the same plot. The 
in-plane d x 2_ y 2 orbital (dashed line in Fig. ^) is predom- 
inantly occupied and shows weak doping dependence in 
almost whole presented range of the hole concentration. 
While axially directed d 3z 2_ r 2 orbital (solid line in Fig. 
|2|) is emptied-out linearly with doping, i.e. all the doped 
hole in this doping range resides on d 3z 2_ r 2 orbital. 
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FIG. 2. Occupancy of d 3z 2_ r 2 (solid line) and d x 2_ y 2 
(dashed line) orbitals as a function of hole concentration. 
The open (filled) circles are the experimental data for d x 2_ y 2 
(d 3z 2_,,2) orbital occupancy from Ref. 4. 




(0,0) (Jt,7l) 
FIG. 3. Bonding 3z 2 — r 2 and x 2 — y 2 bands (in units 
t — 1) referred to chemical potential for hole doping x — 0.3 
(dashed line) and x = 0.4 (solid line). 



lower d x 2_ y 2 band is only weakly changed. 

This doping dependence of the band structure is ex- 
plicitly shown in Fig. |^, where the bonding d x 2_ y 2 and 
d^ z 2_ r 2 bands are presented for hole concentration x — 
0.3 (dashed line) and x — 0.4 (solid line). 

Let us know discuss the model parameters used to re- 
produce the experimental data. The value of the hop- 
ping integral t = 0.18 eV is in the range estimated for 
manganese oxides t ~ 0.1 ~ 0.3 cV.E3 We also point 
out, that some of the authors-estimated hopping inte- 
gral to be higher t ~ 0.7eV.L3 Moderate value of the 
effective inter-orbital Coulomb repulsion U c s = 0.7 eV 
is sufficient to give reasonable fit to the,-data, which jus- 
tifies the MF treatment adopted here.B The AFM su- 
perexchange energy of ti g spins remarkably coincides 
with one (2SJ = lOmcV) £atimated from the spin-wave 
data of Ndo.45Sro.55MnO3.Ll Another important parame- 
ter is the JT splitting of the orbital states, that is dif- 
ficult to directly detected experimentally. For x = 0.3, 
with the above estimate_pf the electron-phonon coupling 
constant (g = 0.5 eV)t3 one obtained the JT splitting 
2Ajt — 0.035 eV. Hence, in the doping regime consid- 
ered here, the JT binding energy is much smaller than 
the carriers band- width. This explains why the pola- 
ronic effects, not considered in the present paper, can be 
ignored. 
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The above discussed evaluation of the orbital occu- 
pancy upon doping can be understood as follows. In the 
doping range x < x c , (x c ~ 0.26 being the hole concen- 
tration at which both orbitals are equally occupied [see 
Fig. |^]) the JT distortion is large and stabilizes the axially 
directed orbital. With doping JT induced gap weakens 
and at x = x c the kinetic energy term supporting the 
in-plane orbital state takes over. Therefore, Sn changes 
sign at x — x c and the d x 2_ y 2 orbital starts to be pre- 
dominantly populated. For x > x c correlation induced 
splitting 2A cor = USn is negative and supports in-plane 
orbital. With farther doping, the orbital anisotropy Sn 
increases, leading to increase of correlation induced split- 
ting. Therefore, the axial band is pushed up relative to 
the in-plane band upon doping, while the chemical po- 
tential is practically pinned. Therefore, all the doped 
holes go to the upper d 3z 2_ r 2 band and its occupancy 
linearly decrease upon doping, while the population of 
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